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This paper studies a scheme of two spatially distant oscillator systems that are connected by Gaussian fields 
and examines distributed entanglement generation between two continuous-mode output Gaussian fields that 
are radiated by the oscillators. It is demonstrated that using measurement-feedback control while a non-local 
effective entangling operation is on can help to enhance the Einstein-Podolski-Rosen (EPR)-like entanglement 
between the output fields. The effect of propagation delays and losses in the fields interconnecting the two oscil- 
lators, and the effect of other losses in the system, are also considered. In particular, for a range of time delays 
the measurement feedback controller is able to maintain stability of the closed-loop system and the entangle- 
ment enhancement, but the achievable enhancement is only over a smaller bandwidth that is commensurate with 
the length of the time delays. 
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I. INTRODUCTION 

Entanglement between quantum systems is considered to 
be an important resource to be exploited in many quantum- 
based technologies proposed in recent years. In particular, 
effective entanglement distribution in a quantum network is a 
problem that has attracted attention in the literature due to its 
importance in applications |1, 2|. However, entanglement is 
very fragile and the amount of entanglement can quickly be 
lost due to decoherence. One strategy to overcome this is to 
use several copies of quantum systems with a limited degree 
of entanglement and to process them to obtain a single copy 
that contains a higher degree of entanglement. This process is 
known as entanglement distillation |3 |. 

Researchers have considered entanglement distillation in 
both discrete and continuous variables. In the continuous vari- 
able case, a particular class of systems of interest are oscilla- 
tor systems that are in a Gaussian state |4|. If one has several 
copies of bipartite entangled pairs of Gaussian states then a 
no-go result of 1 5 1 states that it is not possible to distill further 
entanglement using only Gaussian local operation and clas- 
sical communication (LOCC). Essentially, because the oper- 
ations are LOCC there can be no entangling operations be- 
tween any of the oscillators as they would necessarily have 
to be non-local. However, often in practice one considers dy- 
namical quantum systems having an entangling interaction on 
during the evolution. The effectiveness of entanglement gen- 
eration may be hampered by the decoherence taking place, 
possibly causing the entanglement to eventually vanish. Since 
there is a non-local entangling operation in effect, the no-go 
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theorem does not hold, and this opens possibility to use Gaus- 
sian LOCC operations to protect entanglement during the sys- 
tem evolution. Thus, Mancini and Wiseman |6| considered 
the use of measurement feedback to improve entanglement 
between two cavity modes ai and a2 coupled via a two mode 
squeezing Hamiltonian H = ^(ea^aj — e*aia2); in this paper 
the notation * denotes the the adjoint of a Hilbert space opera- 
tor or the conjugate transpose of a matrix of complex numbers 
or operators. At the same time, the two modes undergo deco- 
herence, due to the cavity photons escaping through the trans- 
missive cavity mirrors on which measurements can be made. 
They showed that there exists a Gaussian LOCC strategy re- 
alized by measurement-feedback that can help to improve the 
amount of entanglement between the two cavity modes com- 
pared to when no measurement and feedback is used. It was 
subsequently shown in | 7 1 that the strategies proposed in |6 | 
are optimal for that particular physical setup. Several other 
works in the literature have also considered using feedback 
for entanglement control, e.g., ISUTTl. 

The main contribution of this paper is developing a 
feedback-controlled scheme that uses distributed resources 
(parametric amplification at two spatially separate sites) 
to generate entanglement between two spatially separated 
continous-mode Gaussian fields. The system of interest is a 
quantum network consisting of two spatially separated open 
Gaussian oscillator systems. They are interconnected via trav- 
elling Gaussian quantum fields that act as common baths be- 
tween the two oscillators and as a source of effective interac- 
tion between them. The scheme also exploits measurement- 
feedback control to enhance the entanglement. It is distinct 
from a measurement-feedback control based scheme on a lin- 
ear quantum system studied in the earlier work |6| in several 
ways: 

1. The scheme uses distributed resources. That is, entan- 
glement is generated exploiting contributed resources 
at two spatially separated locations (say, at Alice and 



Bob's). This use of spatially distributed resources for 
entanglement generation is a key feature of the pro- 
posed scheme. In contrast, | 6 1 considers two oscillators 
interacting in a x*^^^ nonlinear crystal at a single loca- 
tion, and the entangling process occurs only at one site 
using only resources available at that site. 

2. The paper is concerned with entanglement between two 
output fields, that consist of a continuum of modes, 
rather than entanglement between single mode internal 
oscillator modes as in 1 6 1 . Although the system in (61 
does have output fields, they are only provided as ob- 
servables to be measured and fedback rather than left 
freely as entangled resources. In contrast, the present 
work utilizes several output fields, some of which are 
for measurement and feedback while others remain free 
to be used as entangled resources. The rationale for 
considering entanglement in the output fields, rather 
than between internal oscillator modes, is that the out- 
put fields are more easily accessible to be exploited for 
various purposes. For instance, if multiple copies of 
the system are available then the multiple outputs of 
the various copies could, say, be passed onto some (at 
this stage hypothetical) entanglement distiller to pro- 
duce another signal with improved EPR-like qualities. 

3. Since generation of entanglement is via transmission 
channels linking the two spatially separated sites, there 
are inherent features of this scheme that need to be con- 
sidered and which are not present in [6] . The first is the 
presence of losses along these transmission channels, 
and the second is the presence of time delays required 
for the interconnecting fields to propagate between the 
two sites. 

In this work we show that a linear quadratic Gaussian 
(LQG) measurement-feedback controller (Tol [121 [IS can en- 
hance the entanglement between the output fields as compared 
to when there is no controller present. Moreover, this con- 
troller can provide an enhancement even in the presence of 
losses and delays in the overall system. That is, the controller 
displays robustness with respect to the presence of these im- 
perfections. Simulation results are presented and discussed to 
compare the performance of the uncontrolled and controlled 
systems in the ideal case and in the presence of imperfections. 



II. PRELIMINARIES 

A. Linear quantum systems 

In general, the dynamics of an open quantum system, which 
does not contain a scattering process, can be characterized by 
the system-environment coupling Hamiltonian 

m 

i^int(t)=Z^(L,OW*-i^*OW), (1) 

where ^^(t) (j = l,...,m) is the field operator describ- 
ing the j-th environment field and Lj is the system operator 



corresponding to the j-th coupling |[T4l . When the Markov 
limit is imposed on the environment, the field operators sat- 



5{t — s)5i 



A linear quantum system with n-bosonic modes aj {t) {j = 
1, . . . , n) satisfying [a^, a*] = Sij appears when Lj is linear 
and H is quadratic in aj and a*. In this case, the Heisen- 
berg equation of aj{t) = U{tyajU{t) with unitary U{t) = 
ex{) {—i Jq Hi^t{s)ds) has the following form: 



z{t)=Az{t)^B^{t), 



(2) 



where we have defined 



with quadratures = + a*, = {aj — a*)/i, = + 
and = (^j — The field operator also changes to 

^out,j{t) = U {ty S,j {t)U (t) , and measuring this output field 
generates the classical signal 



y{t) = Cz{t)^Dm- 



(3) 



The system matrices A^B^C, and D have specific structures 
to satisfy [a^(t), = Sij for all t >0, and [y{t)^ y{s)] = 

Ofor all5,t > 0. 



B. Measurement-feedback LQG control 

The field ^(t) can closely be approximated by a coherent 
light field generated from a laser device. Thus feedback con- 
trol can be realized by replacing (t) by a modulated field 
Uj{t) + ^j(t), where Uj{t) is a function of y{s) {s <t). 

In the LQG feedback control scheme, the control input 
u{t) = {ul, ixf , . . . , ^1^, u^^) is generated as an output of the 
following classical linear system with y{t) as the input: 



Zc{t) = AcZc(t) + Bcy(t), u(t) = CcZc(t). 



(4) 



Zc{t) is a vector of real c-numbers representing the state of 
controller. Be, and Cc are real matrices to be designed. 
Combining Eqs. ([2]), ([3j, and fij, we have a closed-loop dy- 
namics with variable z = , zj)^. For this system, we 
consider the following quadratic-type cost function: 



^ \^z{tfPz{t) + u{tfQu{t)^dt^ , 



J(u) = lim -E 

(5) 

where P > and Q > are weighting matrices that should be 
chosen appropriately. The expectation is defined as E(X) = 
Tr[X(p(g) |0) (0|)], where p is the initial system Gaussian state 
and |0) is the field vacuum state. The optimal LQG feedback 
control input is given as a minimizer of the cost function ([5]); 
the minimization problem min J{u) can be efficiently solved 
using the Matlab Control System Toolbox to yield the optimal 
set of matrices Ac^Bc, and Cc- 



C. Frequency domain entanglement criteria for 
continuous-mode output fields 

Let us here particularly focus on two output fields ^out,i{t) 
and iout,2{i)' These are continuous-mode fields, therefore 
we need to move to the frequency domain to evaluate their 



entanglement; for the quadratures Xj = ^, 



out,j 



i 



out J 



and 



Vj — i^out,j — Cnt,j)/^' define their Fourier transforms by 
Oj {iuj) = J Oj {t)e-''^^dt/V2^ with o = x, and O = X, Y, 
respectively. Then the two output fields S,out,i{t) and ^ont,2 (t) 
are entangled for the mode at frequency uj if [16J 

V^{iuj)^V-{iuj) <4, (6) 

where V-^{iuj) and V_ {iuj) are defined by the identities 

{{Xi{iuj) ^ X2{iuj)y{Xi{iuj') ^ X2{iuj'))) 

= {{Xi{-iuj) ^ X2{-iu)){Xi{iuj') ^ X2{iuj'))) 

= V-^{iuj)S{uj — uo')^ 
{{Y^iiuj) - Y2{iuj)y{Yi{iuj') - Y2{iuj'))) 

= ((yi(-zc^) - Y2{-i^)){Yi{icu') - FsM))) 

= V-{iuj)5{Lj - co'). 

That is, V-^{iuj) and V-{iuj) are power spectral densities of the 
fields. The inequality ([6]) is a well-known sufficient condition 
for entanglement in the frequency domain [15 , 16|. In the 
case of a two mode-squeezed state, the two power spectra are 
identical and, in an ideal limit, they converge to zero for all 
cj, implying that the so-called Einstein-Podolski-Rosen (EPR) 
pair of fields is produced 1151 . 

III. THE SYSTEM MODEL 
A. Description 



the system Hamiltonian H2 = —e{a2b2 + a2b2)/2\^ with 
the same e > 0. Also a2 is coupled to ^in,2i and ^in,22 via 
L21 = V^^2 and L22 = \/hi/2a2, while 62 is coupled to 
6n,23 via L23 = v^^2- As in Gi, possible losses in the two- 
mode squeezing process in H2 is modeled by coupling a2 and 
62 to additional noise fields (,ioss, 21 and (,ioss, 22, respectively, 
via L24 = ^Jx^2a2 and L25 = ^/x/2h2 with the same con- 
stant X as before. The input fields iin.ji, S,in,j3^ S,ioss,ji^ and 
iioss,j2 for j = 1,2 are in the vacuum state. 

We allow the possibility of photon losses in the two trans- 
mission channels connecting Gi and G2. The losses are mod- 
eled by inserting in the two transmission paths a beam splitter 
with transmissivity a and reflectivity /3, with + = 1. 
Each beam splitter BSj has two ports, one for the incoming 
signal S,out,j3 and an unused port for a noise field ^bsj • More- 
over, we also allow the possibility of time delays along these 
transmission channels that are represented in the figure by the 
(Jt blocks, with T a positive number indicating the transmis- 
sion delay. The time delay block acts on a signal ^ coming 
into the block as ((JtO (^) = ^ ~ Thus the interconnect- 
ing fields satisfy a{aT^out,i3){t) ^ P£,bs,i{'^) = ^in,22{t) and 
<^(crT^ont,23)(^)+/^^B5,2(^) = 6n,i2(^). The dynamics of the 
whole network is then given as follows: 

aiW = -(i + i + f + ^^iW* - V^^-aiW 

X, 



kit) = ^Mty - (y + f - ATCm.iaW 



12(^), 



We consider the system shown in Fig. [T] The spatially 
separated open quantum systems Gi and G2 are connected 
by Gaussian quantum fields. The system Gj consists of two 
oscillator modes aj and bj having the same oscillation fre- 
quency coupled to three independent quantum white noise 
fields (,in,ji, 6n,j2, and ^injs; they are continuous-mode 
and satisfy [£,in,ji{t),£,inji{sy] = ^{t — s). The oscillator 

- ^jk^ 
2. 



modes satisfy the usual commutation relations a^] = 

[^j^K] = ^jk^ [c^jM] = 0, and [a^-, 6^] = for j, k = 1 



In system Gi, the modes ai and hi are coupled via the two- 
mode squeezing Hamiltonian Hi = ie{albl — ai6i)/2>/2 
with e constant. Moreover, ai is coupled to ^in,ii and ^in,i2 
via the coupling operators Ln = y^ai and L12 = Kl2a\, 
respectively, for some coupling constants 7 and k. Also, 
h\ is coupled to ^in,i3 via L\3 = ^/K^bl for coupling co- 
efficient K,i. We allow the possibility of losses in the two- 
mode squeezing process, which is modeled by the interac- 
tion of ai with the additional quantum noise field (,ioss, 11 via 
Li4 = ^Jx^2a\ with x a coupling constant. Similarly, h\ 
interacts with iioss,i2 via L15 = ^J~x/2hl. The system G2 
has a similar structure. The modes a2 and 62 are coupled via 



«2(t) = -{1 + 1 + \)a,{t) + ^h^ity - VliinMt) 

m} 

hit) = ^«2W* - (y + ^)b2{t) - V^UMt) 



^6oss, 22(^)7 

with outputs 



^out,i2{^) = ^CLi{t) + a^/i^{aTb2){t) + a{aT^in,23){t) 

+ /5^B5,2(t), 




(a) 




^lossM - 



-> ^outM Beam splitter 1 ^ 

^out,\2 . , 



Or 



^BSl 

Beam splitter 2 > ^ 
\ 



(b) 

FIG. 1: (a) Input and output fields of the open oscillators systems Gi and G2. In both systems, the output fields corresponding to the loss 
field ^ioss,ij cannot be essentially observed, hence they are not shown, (b) The connected open oscillator systems of Gi and G2- The traveling 
fields contain possible losses modeled by two beam splitters BSl and BS2 with vacuum noises ^bs,i and ^bs,2 entering into their unused 
ports, respectively, ar denotes the operation bringing a time-delay T on the traveling fields. 



^ont, 22 (^) 



a2{t) + aV^(crT^l)(t) + Q^(crT6n,13)(^) 



two spatially separated systems effectively interact through 
this two-mode squeezing Hamiltonian. 



The oscillator systems Gi and G2 could, in principle, be 
physically realized using optical cavities with x*^^^ nonlinear 
crystals, as depicted in Fig. [2] More precisely, the bow-tie 
type cavity Gi contains two modes ai and hi that overlap and 
interact in a x*^^^ nonlinear crystal driven by a classical pump 
beam of effective amplitude e/>/2. The modes ai and hi are 
frequency degenerate but orthogonally polarized. The mirrors 
composing the cavity are partially transmissive, depending on 
polarization of light fields; in particular, the mirrors Mn and 
Ml 2 are partially transmissive for ai but perfectly reflective 
for hi, while M13 is partially transmissive for hi but perfectly 
reflective for ai. The transmittance of Mn, Mi 2, and Mi 3 
are Tn = ^l/c, T12 = i<il/2c, and T13 = kiI/c, respectively, 
with / the optical path length of the cavity and c the speed of 
light. In addition, the field ^in,i2 entering through mirror Mi 2 
and the output iout, 12 are both passed through a 180° phase 
shifter. The system G2 is similarly realized by a bow-tie type 
cavity. For a more detailed description, see e.g., fTT|. 

The particular model described above is of interest because, 
in the large limit of the parameters that adiabatically eliminate 
hi and 62 fTSl, and taking the limit of zero time delays along 
the transmission channels |T9] [201 , the remaining modes ai 
and a2 couple via the Hamiltonian H = —iahi{a\ 
this is a two-mode squeezing Hamiltonian which is the ba- 
sis for generating entangled photon pairs in a nondegenerate 
optical parametric amplifier (NOPA) 1 15 |. In this sense, our 
system is a realistic approximation to the ideal system where 



B. Quadrature form and transfer function 



We now specialize to the case where the transmission delay 
T is negligible compared to the time scale of the dynamics 
of the systems Gi and G2. Then, in terms of quadratures, 
the coupled network Langevin equations with no transmission 
delays are: 



I 2^in,23 - 2^BS,2 



KKi 



in, 11 



Zoss, 11' 



BS,2 



2 



^Zoss, 11' 



2 ^loss,12i 



2^/2 



^1 ~ (y + f )^1 ~ V^C>i,13 - V f ^Ls,12' 




^ out, 11 



FIG. 2: Candidate for a realization of the systems Gi and G2 in quantum optics. The modes ai and hi are frequency degenerate but 
orthogonally polarized. The mirrors composing the cavity Gi are partially transmissive, depending on polarization of the light fields; the 
mirrors Mn and Mi 2 are partially transmissive for ai but perfectly reflective for hi, while Mi 3 is partially transmissive for hi but perfectly 
reflective for ai. The cavity G2 is also constructed in a similar way. For simplicity, other optical devices such as phase shifters are not shown. 
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hi 
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2 ^/oss, 21' 



-<i + i + i>< 



13 
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^Zoss, 21' 



6« = 

2V2 



(y + f - ^/^Cl 



in, 23 



^l^fn,23 



SZoss, 22' 



X^p 

2 ^^oss, 22' 



with outputs 



out jll 



Q _ 



^^ont, 12 
tP 

Sont, 12 



Sont, 21 



Sont,22 



out,22 



V7«? + C'^n,ll' Cn^,ll = V7«? + Cll' 
-^1^1 + «V^^2 + <n,23 + /^^l5,2' 

+ C'^n,21' Cn^,21 = ^7^2 + ^,21' 



Now, let us observe some properties of the above dynamical 
equations for the oscillator and field quadratures. In particular, 
from the above equations it can be verified that the equations 
are not fully coupled and in fact: (a) The set of equations 
for al^a^jbl, and form a closed set of equations driven 
by the commuting set of noises ^f^ 23' ^i^n,i3' ^?n,2i' 

^Ls,ll' iloss,12^ iloss,2V ^Ls,22' ^B5,l' ^B5,2' 

the set of equations for a^^a^jb^, and b^ form another closed 
set of equations driven by the commuting set of noises ^f^ n, 

dP tP tP tP tP t^l tP ' 

^in,23' ^in,13' ^in,21' ^^oss, 11' ^^oss, 12' ^Zoss, 21' ^'loss, 22' SBS',1' 

and 2- Although we are considering the case of no time 



delays, it may be easily inspected that the decoupling between 
the above sets of closed equations for certain quadratures and 
commuting noises also holds in the time delay case since the 
structure of the equations are precisely the same. 

We can now consider the Fourier transforms of the observ- 
ables above and the system transfer functions. Introduce the 
column vector of system operators 

z,{t) = {al{t),al{t),bHt),bP{t),)^, 
Z2{t) = {a\{t),aP^{t)y,{t),bl{t)f, 

and 



_ {dQ tP dQ dQ 
1 ~ VSin, 11' Sin, 23' Sin, 13' Sin, 21' 



dQ dQ dQ dP dQ dP \T 

SZoss, 11' SZoss, 12' SZoss, 21' SZoss, 22' ^55,1' ^55,2/ ' 

d — (dP dQ. dP dP 

S2 — VSin, 11' Sin, 23' Sin, 13' Sin, 21' 

dP dP dP d^ dP d^ \T 

SZoss,ll' SZoss,12' SZoss,21' SZoss,22' SBS',1' SBS',2/ ' 

?=(Ci(tf,6(irr. 

By the decoupling structure already noted above, we have 

Zi{t) = A,Zj{t)-rBjij(fy, (i = l,2), 

Ct,ii(*) - Ct,2i(o = ^^2^2(t) + D^m. 

where Ai^A2jBi, and B2 are real matrices of suitable di- 
mensions whose entries can be readily determined from the 
equations for aj, a^, and b^, and the row vectors Ci, C2, 
Di, and D2 are given by: 

Ci = (l, 1, 01x2), C2 = (l, -1, 01x2), 

Di = (1, 0ix3, 1, Oixe), D2 = (1, 0ix3, - 1, Oixe). 

Let xj{t) = CutA^) and Vjit) = CutA^) and let X„ 
Yj, and denote the Fourier transforms of Xj, yj, and 
respectively (see Section [nc]). Then 



1,2). 



Xi{iuj) -\- X2{ioj) = Hi{iuj)Ei{iuj)^ 
Yi{iuj) - Y2{iuj) = H2{iuj)E2{iuj), 



with Hj{iuj) = Cj{iujl 

Using the fact that {{Ej{iujy)^Ej{iuj')^) 



(Ej{—iuo)'E.j{iuo')^) = /ioxio^(<^ — ^0' we find 
that {{Xi{iuj) + X2{iuj)Y{Xi{iuj') + X2{iuj'))) = 
TT[Hi{iujyHi{iuj)]S{uj - uj') and {{Yi{iuj) - 
Y2{iuj)y{Yi{iuj') - Y2{iuj'))) = TT[H2{iujyH2{iuj)]S{uj - 
uj'). Thus, for the uncontrolled network we conclude that 

V-{iuj) = TT[H2{iujyH2{iuj)]. 



IV. LQG FEEDBACK CONTROL OF THE SYSTEM 

The control scheme is shown in Fig. [3] The input to the 
controller will be the signals yc,ii, ^c,i2, ^c,2i, and ?/c,22 that 
are obtained by performing dual homodyne detection on the 
output fields ^ont, 12 and (,out,22, respectively. Two dual ho- 
modyne detectors, labeled 1 and 2, are required, and each 
consists of a 50:50 beam splitter, where at the unused beam 
splitter port a vacuum noise source ^^^j, {j = 1, 2) comes in. 
At the output of one beam splitter, the position quadrature of 
the field is measured while at the other output the momentum 
quadrature is measured. For the j-th dual homodyne detector, 
the outputs i/cji and ycj2 are given by: 



a 



V2 



(7) 



(8) 



1 



V2 



(9) 



are u^^^it) + Cl^i and vF^^(t) + Cf„ji, respectively, while 
the position and momentum quadratures of the output field 
M0Dj2 are u%(t) + ^Lj-s and u]^{t) + ef„,,3' respectively. 
Lett/c(i) = (yc,ii(i),yc,i2(i),yc,2i(i),yc,22(i))^ and 

/ q p q p q p q p \T 

'^c — I'^li 7 '^117 '^21' '^21' ^12' ^12' '^22' '^22/ • 

The LQG controller has an internal 8-th order state Zc{t) that 
obeys the following classical Langevin equation: 

Zc{t) = AcZc{t) + BcVcit), Uc{t) = CcZc{t), (11) 

where Cc are real matrices of the appropriate dimen- 

sions. Here, the order of Zc{t) is 8, since the degree of the 
system to be controlled is 8, corresponding to the oscillator 
quadratures aj, a^, and 6^ for j = 1, 2. 

We assume that the controller is located on the site of G2 so 
that the delays in transmitting the control signals U21 and U22 
from the controller to the system G2 are negligible. However, 
we allow the possibility of delays in transmitting the control 
signals uu and U12 from the controller to Gi ; here we assume 
those time delays take the same quantities and let us denote 
them by T^. Then, the dynamical equation for the closed loop 
system can be obtained simply by making the substitutions of 
Eqs. ([7])-(p^ into Eq. ( pT) and the substitutions 

O'Tmi'^ll + 6n,ll) ^in,ll, CFt^{ui2 + 6n,13) ^ 6 
U2I + 6n,21 6n,21, '^22 + ^in,23 ^ ^in,23 

into the dynamical equations of the system. Now define 

Z — \Z ^ Z^ ) , — 5 S/i,l7 S/i,2' S/i,2>' • 

The closed-loop system without time delays is then given by 

k{t)=Az{t)^Bi{t), 
where A and B are real matrices of the form: 



an,13 7 



A 



Ai Ai3 
A2 A23 
A31 A32 A33 



B 



5i 

^2 

Bz\ Bs2 Bs3 



a 

1 

71 



«n,23)W 



V2 



(10) 



The controller produces 8 control signals u^i{t), Uji{t), 
1x^2 (^)' and u^2i'^) J ~ ^ which drive the quantum sys- 
tem via modulators. The control signals u^^ and u^^ form 
the real and imaginary component of the complex classical 
signal Uji — u^^ + iu^^ that will drive the modulator la- 
beled M0Dj7 for j, / = 1,2. The output of the modulator 
M0Dj7 then drives the system Gj, see Fig. [3] The position 
and momentum quadratures of the output field of MODjl 



Here, Ai, A2, Bi, and B2 are as defined previously in Sec. 



IIIB and Aji and Bji are matrices that are determined by 
the resulting closed-loop system. The equations for ^^^^ + 

Cut,2i and ^^^^11 - ^^^^ 21 are now of the form: 

Cnt,llW +Ct,2lW = (Cu Oix4, Ci2)^(t) + AlW, 
CtaiW -C21W = (Olx4, C2, C22)z{t)^D2m, 



with Ci and C2 also as given in Sec. |IIIB 
{Di, O1X14) and£>2 = (Oixio, ^2, Oix4 



and Di = 



Note that there 

are contributions from the classical controller state Zc{t) of 
the controller in these quadratures. However, since we are in- 
terested in the entanglement between the output fields, we can 
omit these contributions, since they, being classical, have no 
bearing on the degree of entanglement. Thus, in the case of 




FIG. 3: The controlled quantum network. The LQG controller is located on the site of G2; hence the time delays of communication between 
the controller and the system G2 is negligible, while we assume that the control signals u\\ and u\2 propagate from the controller to Gi with 
time delays Tm- The operation bringing this time-delay is denoted by ctt^ . 



a measurement-feedback controller in the loop it suffices to 
consider the modified outputs 

Ct,llW -Ct,2lW = (Olx4, C2, 0ix8)^W +^2^». 

Define = C,,nW' ^^{t) = C^^M^), yi{t) = 

CutM^)^ and y2{t) = C^^M^) and let X,{iuj), X^iiuj), 
yi(ico'), and Y2{iuj) be their Fourier transforms, respectively. 
Also, let E{iuj) denote the Fourier transform of ^{t). Then, 
the transfer function Hi {iuj) from 'E{iuj) to Xi {iuj) + X2 {iuj) 
and H2{i(jj) from to Yi{i(jj) — Y2{i(jj) are given by 

Hi{iuj) = (Ci, 0ixi2)(^c^/ - A)-^B + ^1, 
H2{iuj) = (Oix4, C2, 0ix8)(^c^/ - A)-'B + D2. 

Then, noting that {{E{iujy)^E{iuj')^) = 
{E{—iuj)E{iuj')^) = (/24x24 + iZ)5{uj — uo'), we ob- 
tain the expressions: 

((Xi(ic^) +X2M)*(Xi(icjO +^2(^C^0)) 

= Tr[(^i(zc^)*^i(zc^))^(/24x24 + - 



((FiM - Y2{iuj)r{Yi{iuj') - Y2{iuj'))) 

= Tr[(^2M*^2M)^(/24x24 + iZ)]S{uj - UJ'), 
= Tr[^2M*^2M]^^(^ - ^0- 

Here, Z is a real invertible 24 x 24 skew-synmietric ma- 
trix given by Z = — icj)S(i(x;)^)^). 
Therefore, we conclude that the expressions for {iuj) and 
V- when there is an LQG controller in the loop are 

V+{iuj) = TT[Hi{iujyHi{iu)], 

V-{iu) = Tr[^2M*^2M]. 



V. ENTANGLEMENT BETWEEN CONTINUOUS-MODE 
GAUSSIAN OUTPUT FIELDS: THE IDEAL CASE 

We here consider the entanglement generated between 
Cout, 11 (^) and £,out,2i{t), with and without an LQG feedback 
controller, in the idealized situation where there are no losses 
in the two mode squeezing processes (i.e., x = 0) there 
are also no losses in the transmission channels between Gi 
and G2 (i.e., a = 1). First, we assume that the transmission 
delays along the transmission channels are negligible. Then 
we will consider how the control performance and the closed- 
loop stability are affected by non-negligible time delays. 



A. Negligible transmission delays 

Throughout we will consider the case where k, = 1.8 x 10^ 
Hz, 7 = 1.5 X = 2.7 X 10^ Hz, /^i = 10 x = 1.8 x 10^ 
Hz, and e/V2 = ^/kk^ = 4.0249 x 10^ Hz. In the case of 
the optical system shown in Fig. [2j these parameter values are 
realized using the mirrors with transmittance Tn = T21 = 
0.045, T12 = T22 = 0.015, and T13 = T23 = 0.3, where the 
optical path lengths of each cavities are both set to / = 0.5 m. 
In this section we also set a = 1 and T = Tm = 0. To design 
an 8th-order LQG controller, let us set the cost function ([5j in 
the following form: 



J{Uc) 



lim -E 

T^oo T 



^ {g{[Cr C2]z{t))' 
Uc{tfuc{t)}dt^. (12) 



The weighting constant g is taken to be ^ = 1 x 10^. With 
the use of the Matlab Control System Toolbox, we obtain the 
optimal LQG controller, and we show the frequency domain 
power spectra plots V-^{iuj),V-{iuj), and V-^{i(jo) -\-V-{i(jo) in 
Fig.|4]in dB scale (here A > in linear scale is 10 logio ^ 
dB scale). These figures indicate that all three power spec- 
tra are lower when the LQG controller is present as com- 
pared to when there is no controller. Here the controller 
can provide an additional attenuation to all three spectra by 
slightly more than 2.6 dB up to frequency of about 10^ 
rad/s. Moreover, by the entanglement criterion ([6]), we see 
that entanglement is achieved for modes with frequency up 
to slightly above 10^ rad/s but less than 10^ rad/s (note that 
entanglement is achieved at a mode of frequency u whenever 
V^{iu) + V-{iuj) < 101ogio(4)dB = 6.0206 dB). Note that 
this is achieved despite the controller being designed to mini- 
mize the cost function ([12]) rather than minimizing any of the 
three power spectra directly, indicating the utility of the cost 



Slightly above 10^ rad/s, there is no more reduction and in 
fact one can see a marked increase in the magnitude of each 
spectra at certain frequencies above 10^ rad/s. This suggests 
that the feedback controller is still effective for enhancing the 
entanglement between the output fields even in the presence 
of time delays, but the bandwidth at which enhancement is 
achieved is reduced. However, additional care has to be taken 
before we can be conclusive about this claim. As we have 
seen, there are frequency ranges in which all three spectra 
experience a large increase, which could be indicative of in- 
stability (instability here is in the sense that the symmetrized 
covariance matrix of the system diverges as t 00). Lin- 
ear delay differential systems can be described by abstract 
infinite-dimensional differential equations and numeri- 
cal algorithms are available to examine the stability of these 
systems. Here we use the freely available DDE-BIFTOOL 
toolbox (2211231, a Matlab toolbox to determine the stability 
of a delay differential system, and find that the system under 
consideration is indeed stable, the right hand most root of the 
characteristic equation of the system has a negative real part 
|29|. 

Note that the log-log nature of the plots and the number of 
points used to produce Fig. [5] (as well as Fig. 10) in Matlab 
give the impression of non-smooth functions, but in fact this 
is not the case. Since the analysis shows that the closed-loop 
system is stable, the closed-loop transfer function has no poles 
on the right half plane and all poles are bounded away from 
the imaginary axis and thus the functions V±{iuj) depicted in 
the plots are theoretically guaranteed to be smooth functions 
of (i.e., they infinitely differentiable functions of cj). How- 
ever, the plots indicate that when there are time delays these 
functions fluctuate faster at some high frequencies. The fluc- 
tuations in the figures are prominent since the time delays have 
been taken to be comparable to the time scale of the node dy- 
namics. If the delays are gradually decreased to zero then the 
high frequency fluctuations gradually smooths out. 



function ( 12 ) for feedback control of output entanglement. 



B. Non-negligible transmission delays 



VI. EFFECT OF AMPLIFICATION LOSSES AND 
TRANSMISSION LOSSES 



Since the dynamics of z{t) and Zc{t) preserves the Gaussian 
nature of the state, we can study an equivalent linear classical 
delayed differential Gaussian system with finite time delays 
(in the sense that the mean and (symmetrized) covariance ma- 
trix of the quantum system and its classical equivalent evolve 
in an identical manner). See, e.g., the text |21 1 for a treat- 
ment of classical linear systems with time delays. Such linear 
differential systems with delays can be handled with the Mat- 
lab Control System Toolbox's ' delay ss' object. We consider 
the case where T = 10"^ s and = 2 x 10"^ s, which 
is about an order of magnitude longer than the time scale of 
the system dynamics. Using the controller that had been de- 
signed in Sec. V A we obtain plots of V-^{iuj), V-{i(jo), and 
V-^{iuj) + V-(iuj) as shown in Fig.|5] They indicate that de- 
spite the time delays there is still a significant reduction of 
about 2.6 dB in all three spectra for frequencies up to about 
10^ rad/s and 1 dB for frequencies between 10^ and 10^ rad/s. 



Here we consider the more realistic case where there are 
losses in the two-mode squeezing process and along the trans- 
mission lines connecting Gi and G2. We will design our LQG 
controller based on the assumption that we do not how much 
these losses appear in the system, effectively. Thus we set 
the controller to be identical to the one we had previously de- 
signed under the assumption x = and a = 1. 

We consider first the case where the amplification loss co- 
efficient is X = 1.3975 X 10^ Hz and there are no transmission 
losses (a = 1). The results are shown in Fig.|6] It can be seen 
that in the frequency region where the controller can reduce 
the power spectra, the reduction is smaller than if there were 
no amplification losses (about 1.4 dB reduction compared to 
about 2.6 dB reduction in the latter). For the same amplifica- 
tion loss coefficient, the cases with transmission losses of 3% 
(a = 0.97) and 5% (a = 0.95) are shown in Figs, [t] and [s] 
respectively. Comparing Fig. |4] and Fig.[6j we see that, as can 
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FIG. 4: Plots of V+ (iu) (left), (icj) (middle), and V- (iu) + V+ (icj) (right) in dB against the frequency in rad/s, without an LQG controller 
(solid line) and with an LQG controller (dashed line) designed to minimize the cost function \12\ . The same lines are used as legends in 
subsequent figures as well. 
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FIG. 5: Plots of V+ (iuj) (left), V- (iuj) (middle), and V- (iuj) + V+ (iuj) (right) in dB against the frequency in rad/s when there are time delays 
present (T = 1 x 10"^ s and = 2 x 10"^ s). 
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FIG. 6: Performance of the system when the amplification loss is 
X = 1.3975 X 10^ Hz, no transmission losses, and there are no 
delays. 



a = 1. That is, the controller exhibits a level of robustness in 
its ability to improve the system performance. 
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be expected, the presence of amplification loss has an adverse 
effect on the EPR-like entanglement that can be observed at 
the output; that is, all the power spectra are amplified at fre- 
quencies up to 10^ rad/s. Figs.[6][8]then show that the presence 
of increasing transmission losses leads to a corresponding in- 
crease in the power spectra across the same frequencies and 
thus a worsening of the quality of the EPR-like entanglement 
in the output fields. However, in all of these cases, it is clear 
that the presence of a controller leads to an improvement in 
the entanglement. Moreover, this is despite the fact that the 
controller was designed under the assumption that X = and 



FIG. 7: Performance of the system when the amplification loss coef- 
ficient is X = 1.3975 X 10^ Hz, transmission losses of 3%, and there 
are no delays. 

Finally, Fig. |9] shows the case where x = 5.5902 x 10^ 
Hz, which is a substantial percentage of amplification loss 
compared to pump intensity, and the transmission loss is 5%. 
Here, again, the time delays are assumed to be zero. The fig- 
ure indicates that the LQG controller still manages to improve 
entanglement although the power spectra are, as can be ex- 
pected, noticeably higher than those in Figs. |4][8] Finally, 
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time scales that are shorter than the time scales for propaga- 
tion of fields between nodes. Thus control of quantum net- 
works with time delays is a timely research topic since this 
scenario would be encountered in many quantum networks of 
interest. To the best of the authors' knowledge, only a few pa- 
pers have so far considered quantum feedback control in the 
presence of time delays, e.g., I 27n28l . 



FIG. 8: Performance of the system when the amplification loss coef- 
ficient is X = 1.3975 X 10^ Hz, transmission losses of 5%, and there 
are no delays. 





FIG. 10: Performance of the system when the amplification loss co- 
efficient is X = 5.5902 X 10^ Hz, transmission losses of 5%, and 
there are time delays T = 10~^ s and = 2 x 10~^ s. 



FIG. 9: Performance of the system when the amplification loss coef- 
ficient is X = 5.5902 X 10^ Hz, transmission losses of 5%, and there 
are no time delays. 



VII. CONCLUSION 



for the same values of x and a, the effect of the presence 
of the time delays T = 10"^ s and = 2 x IQ-^ s is 
shown in Fig. 10 As we had seen in Fig. |5] the delays re- 
duce the frequency range in which a reduction in the power 
spectra by the controller can be observed. In particular, in the 
frequency range up to about 10^ rad/s, where the reduction 
is observed, the amount of reduction is roughly the same as 
what can be achieved without the delays. Also, again using 
DDE-BIFTOOL, we can inspect that the system is stable (i.e., 
the symmetrized covariance matrix of the system converges 
as t ^ oo). Also, although not shown here, we remark that 
in general the controlled system is able to remain stable even 
when the time delays T and Tm are increased up to 0. 1 and 0.2 
seconds, respectively. However, of course, it is undesirable to 
work in such severe cases of time delays, since it means that 
the system reaches steady state rather slowly. 

We conclude by remarking that in the non-quantum setting 
there are sophisticated methods for designing controllers be- 
yond the LQG that specifically take into account time delays 
(see, e.g., |26| and the references cited therein ), which have 
the potential to be adapted to the quantum setting. These more 
advanced controllers may potentially offer further improve- 
ments to the power spectral profile but are, however, beyond 
the scope of the present paper. Time delays are important in 
the quantum network setting where nodes can have dynamical 



This paper has developed and studied a distributed en- 
tanglement generation scheme for two continuous-mode 
output Gaussian fields that are radiated by two spatially 
separated Gaussian oscillator systems. It is shown that a 
LQG measurement-feedback controller can be designed to 
enhance the EPR-like entanglement between the two output 
fields across a certain frequency range, even in the presence 
of important practical imperfections in the system. It is 
demonstrated that the controller displays a degree of robust- 
ness in the sense that although it was designed for the ideal 
scenario it can still provide an enhancement and stability over 
the case of no controller being present, despite the presence 
of the imperfections. In summary, the results reported here 
indicate the potential utility of feedback controllers in the 
task of distributed entanglement generation using distributed 
resources. 
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